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Abstract 

Helium was first isolated on Earth in 1895, by Sir William Ramsey. 
One hundred and one years later, it seems like a good time to review 
our current theoretical understanding of the helium atom. Helium 
has played an important role both in the development of quantum 
mechanics and in quantum field theory. The early history of helium is 
sketched. Various aspects of the modern theory are described in some 
detail, including (1) the computation of fine structure to order a 2 Ry 
and a 5 Ry, (2) the decay of metastable states, and (3) Rydberg states 
and long-range forces. A brief survey is made of some of the recent 
work on the quantum field theory of He and He-\ike ions. 
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1 Introduction 



The title of my talk is inspired by that of a famous novel by Gabriel Garcia 
Marquez: "100 years of solitude." Helium, whose existence was not even 
suspected till the middle of the 19th century, has experienced exactly 101 
years of attention since its terrestrial discovery in 1895. It has played a major 
role in the development of atomic, nuclear, and condensed matter physics. 
The helium atom, and its cousins, the if e-like ions, continue to be a subject 
of active study. It has played an important role, second only to hydrogen, 
both in the development of quantum mechanics and that of quantum field 
theory (QFT). Since this atom is, in the sense of the title of this conference, a 
quantum system par excellence, I thought this would be an appropriate time 
to review the history of this subject, with emphasis on the QFT aspects. 

I will begin (Sec. 2) with a historical sketch concerning the discovery of 
helium, first in the Sun and then on Earth, including some fascinating facts 
which I gleaned by browsing through a wonderful book by Isaac Asimov.jl]] 
(I take special pleasure in referring to Asimov, for several reasons: He was 
born in Belarus (Petrovici, 1920) and he is a distant relative of one of the 
members of the organizing committee! Moreover, he obtained his higher de- 
grees at Columbia University, where I was a student many years later.) I will 
then describe the importance of the role played by study of He in the devel- 
opment of quantum mechanics (Sec. 3) and in the early stages of quantum 
electrodynamics (QED) (Sec. 4). I then turn to developments beginning in 
the 1950's, including the first calculation of atomic energy levels to accuracy 
a 3 Ry (Sec. 5). Next, I will consider the theory of the decays of excited states 
of He and ife-like ions, especially of metastable states, in which the Sun 
enters our story once again (Sec. 6). In Sec. 7, I note the role played by 
study of high-lying, so-called Rydberg states (n > 10) of He in advancing 
the theory of long-range forces arising from two-photon exchange. I conclude 
in Sec. 8 with a brief survey of some of the recent work involving QED and 
ife-like systems. 
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2 The discovery of helium: From Janssen to 
Ramsay 



The story begins with someone whom you have probably never heard of, one 
P. J.C. Janssen (1822-1907). From his name you might think he is yet another 
one of those Scandinavians who have contributed so much to the discovery of 
the elements ... but you would be wrong. The initials stand for Pierre Jules 
Cesar and his name is pronounced in the French manner (Zhahn-sen'). Born 
in Paris, Janssen was a great traveler in the interests of astronomy. He went 
to Peru and the Azores, to Italy and Greece, to Japan and Siam. In 1868, 
still without benefit of United/Lufthansa and the like, he traveled to India 
to observe a total eclipse of the sun. More about that in a moment. But to 
give you an idea of Janssen's daring and determination, I add that in 1870, 
when Paris was besieged by the Prussians, and Janssen was eager to get to 
Algeria to observe another total eclipse, he escaped from Paris by balloon 
flight! During his observations in India, Janssen noticed a strange spectral 
line and forwarded his data to Sir Norman Joseph Lockyer (1836-1920), a 
British astronomer who was a pioneer in the study of solar spectra and coined 
the term "chromosphere". Lockyer compared the position of Janssen's line 
to lines of known elements and concluded that it must belong to an as yet 
unknown element, possibly not even existing on earth. He named it helium, 
after the Greek work for the sun, "helios" . But his suggestion was dismissed 
by chemists, who were not prepared to accept the existence of a new element 
on what was then regarded (not without reason) as a flimsy basis. Indeed 
subsequently, many new lines were discovered and incorrectly attributed to 
new elements, such as "coronium", "nebulium," and so on; all of these turned 
out to correspond to ionized states of known elements. 

Our story now turns to the two R's, Rayleigh and Ramsay. Lord Rayleigh 
(1842-1919) is known to most of us for his work in mathematical physics. But 
he was also a great experimentalist. During an extensive series of measure- 
ments on the density of gases, undertaken as a test of Prout's hypothesis, 
he noted that the density of nitrogen obtained from the air appeared to be 
very slightly higher than that of the nitrogen obtained from chemical com- 
pounds, unlike the case with oxygen. In frustration he wrote to the journal 
Nature, (founded in 1869 by Lockyer himself and still going strong) for sug- 
gestions. In response, after getting permission from Rayleigh to do so (!), the 
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matter was taken up by Sir William Ramsay (1852-1916) in 1892. Ramsay 
remembered that a century earlier Cavendish had found that when he tried 
to combine the nitrogen of the air with oxygen a final bubble of air was left 
over, which Cavendish speculated to be a gas heavier than nitrogen which 
did not combine with oxygen. Ramsay repeated the experiment of Cavendish 
and found a similar bubble. Using the now available spectroscope, Ramsay 
and Rayleigh studied the lines from the left-over gas and identified for the 
first time a zero-valence element, which they called "argon," from the Greek 
word for inert. With the periodic table as a guide, Ramsay guessed that 
there might be a whole column of inert elements and began the search. In 
1895, following up earlier work of the American geologist W. Hillebrand, he 
isolated a gas from a uranium mineral called clevite, previously mistaken for 
nitrogen, which showed the same bright spectral line seen by Janssen over 
25 years before. Helium had been found on Earth! Later, starting with 
liquid argon, Ramsay isolated three more inert gases: neon (new), krypton 
(hidden) and xenon (stranger). In 1904, Ramsay and Rayleigh received the 
Nobel Prizes for chemistry and physics, respectively, for their work on the 
inert gases. 

3 Helium and quantum mechanics: From Bohr 
to Hylleraas 

Within a decade or so after its discovery, helium began to be of great impor- 
tance in the development of atomic, nuclear and condensed matter physics. 
With regard to condensed-matter physics, in 1908 Kammerlingh-Onnes liqui- 
fied helium at 4 K; in 1911 he used helium in the discovery of supercon- 
ductivity. With regard to nuclear physics, helium was in at the beginning: 
Rutherford's discovery of the first nucleus, that of gold, involved an a-particle 
beam and in 1906-1909 Geiger and Rutherford identified the a-particle as the 
doubly ionized helium atom. 

The discovery of the atomic nucleus led directly to the Bohr model and 
to Bohr's explanation of the spectrum of hydrogen in terms of quantization 
of the action. However, the attempt by Bohr and Sommerfeld to extend 
the model to account for the spectrum of helium failed. In the 1920's this 
failure led in part to the discovery of quantum mechanics by Heisenberg and 



4 



Schrodinger and, in particular, to the Schrodinger equation for the stationary 
states of the hydrogen atom or hydrogen-like ions: 



Kmi) = [^ + u c {i)W) = w<i>{i). (i) 

with UciX) the Coulomb interaction of the electron with a point nucleus: 
UciX) = —Ze 2 /4:7rri. With the constraint that is normalizable, one recovers 
the Bohr levels as eigenvalues of h(l). 

What about helium? Within the Schrodinger framework the energy levels 
ought to be determined by 

^(1,2)0(1,2) = W0(1,2), (2) 

where 

/i(l,2) = /i(l)+/i(2) + E/c(l,2), (3) 

with £7(7(1,2) = e 2 /47rr 12 . Since (§) can no longer be solved analytically, 
approximation methods must be used. This problem was attacked in the 
later 1920's by several authors but most successfully by Egil Hylleraas @ 
who, many years before Mandelstam, introduced a set of 5, t and u variables 
into physics which were to become famous: 

s = n + r 2 , t = r x — r 2 , u = r 12 . (4) 

For spherically symmetric (L = 0) states, the spatial wave function 0(1, 2) 
may be taken to depend only on these rotationally invariant variables. As 
an ansatz for the application of the variational principle 

S < | h(l,2) | >= 0, (5) 

Hylleraas wrote 

= 4>(s,t,u) = f(ks,kt,ku). (6) 
Here A; is a "nonlinear parameter" and / is of the form [in atomic units] 

f(s,t,u)=e- s / 2 P(s,t,u), (7) 

with P(s,t,u) a multinomial in s, t, and u: 

P{s,t,u) = Y,c nt2l , m s n t 2l u m (8) 
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involving "linear parameters" cy^. With a six-parameter P, his result for the 
energy of the helium ground state gave for the ionization energy, J t h = 1-8065 
Ry, to be compared with the experimental value at that time, J exp = 1.807 
Ry. This agreement, to better than a part in 10 3 , was important in putting 
the approach to atomic spectra based on a many-body Schrodinger equation 
on firm, quantitative ground. 

In the subsequent 60 years, and especially after the advent of modern 
computers, the increasingly precise study of the discrete spectrum of the 
non-relativistic Hamiltonian ([]]) developed into a small industry. As I will 
describe later, in recent years this study has in fact taken a great leap forward. 
The obvious generalization of @ to an N-electron atom or ion (N > 2) , 
when supplemented by the antisymmetry principle and electron spin wave 
functions, constitutes the foundation of the nonrelativistic theory of atomic 
spectra. 



4 Helium and QED: From Darwin to Breit 

4.1 The patchwork era 

Of course, no matter how accurately (§) is solved, the energy levels obtained 
from it are not accurate enough for comparison with very precise experiments. 
Even in classical electrodynamics (CED), the interaction between two point- 
like charged particles is not given just by their Coulomb interaction, which 
is all that is included in @. Already in 1920, Darwin||] had shown how the 
inclusion of relativistic effects in CED leads, to order v 2 /c 2 , to an effective 
Lagrangian interaction term of the form 

(eie 2 /47rc 2 )(vi • v 2 + vi • rv 2 • f) , . 

Unt Yr • (yj 

Upon quantization, this leads to a version of what later came to be called 
the orbit-orbit interaction, U - . Further, after the discovery of electron 
spin and the associated magnetic moment, it was clear that to describe the 
fine-structure accurately one would have to include not only the spin-self- 
orbit interaction U s - s - (i) of each electron with the nucleus, as in the 
case of hydrogen, but also a spin-other-orbit interaction C/ s _ _ (1,2) 
and a spin-spin interaction C/ S _ s (l,2). These terms could be, and were, 
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written down by analogy with the corresponding classical terms, although 
not completely correctly in the case of U s - S . It was also clear that at this 
level of accuracy (a 2 Ry) one should include, again as in H, the relativistic 
kinetic energy correction Ukin proportional to pf/m 3 , for each electron "i". 

Is that all? The discovery of the (one-electron) Dirac equation led one 
to expect additional terms. Reduction to nonrelativistic form of the Dirac 
equation for H or an H-like ion, viz. 

(Ml) + U c (l))ip = E^/j, [Ml) = «i • Pi + Pim] (10) 

leads not only to the spin-orbit and kinetic energy corrections but also to 
a contact term U con t(l), the so-called "S-state interaction", which has no 
classical analogue. Clearly, there must be counterparts of this in the case 
of helium. In summary, as a result of all this patchwork, one expects the 
leading correction to the energy levels obtained from @ to take the form 

=<^sp\uif f \^ S p> (11) 
where ^sp(l,2) is a Schrodinger-Pauli (SP) type wave function and 

V% = U sA . + U s . d , (12) 
with U s .i. and U s ,d. the spin-independent and spin-dependent parts of the 

(2) 

effective interaction U\tf 

U sA . = nU km (i)+U cont (i)} + £/ -o(l,2), (13) 

U s . d . = EC/ s _ (i) + fW (l,2) + C/ S _ s (l,2). (14) 
I will write down the explicit expressions later [Sec. 5.2]. 

4.2 The next stage 

Even if the interaction between two electrons were perfectly described by 
Ucij-12) , it was evident that to take into account relativistic and spin- 
dependent effects in a unified and systematic way it would be desirable to 
construct a generalization of the one-electron Dirac equation to two elec- 
trons. The first steps in this direction and in improving the above patchwork 
approach were taken by J.B. Gaunt J4[ and especially by Gregory Breit.0 
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4.2.1 The two-electron Dirac-Coulomb Hamiltonian. 

As a starting point for a relativistic description of the He atom, Breit 
took what appears to be an obvious extension of (|10D : He replaced the 
4— component Dirac spinor ip{l) by a 4 x 4 = 16 component spinor ip(l,2) 
and required that it be an eigenfunction of an operator which is simply the 
sum of the h'ffii) , and Uc(l, 2) . I like to call this operator the (two-electron) 
Dirac-Coulomb Hamiltonian and denote it by Hdc0->2) '■ 

H DC (1,2) = ht{l) + ht{2) + U c (r 12 ). (15) 

In this notation, Breit 's requirement takes the form 

H DC ip(h2) = £ty(l,2). (16) 

The normalizable solutions of this equation were presumed to give a good 
description of the discrete states of the helium atom. 

4.2.2 Transverse photon exchange and the Breit operator. 

Next, Breit took into account effects arising from the interaction of the elec- 
trons with the quantized radiation field A^(x) by including the interaction 
which, in this context is given by 

H T = Sea* • A T (r<) (17) 

Using second-order perturbation theory and dropping terms which corre- 
spond to self-interaction, Breit argued that to a good approximation the 
level shift AE arising from this interaction would be given by 

AE=<^\U B \ip>/<^\^> (18) 

where Ub is an operator which came to be known as the Breit operator, 

U { b' 2) = -(e 2 /47r)(a r a 2 + a r ra 2 ■ f)/2r, (r = n - r 2 ) (19) 

and ip is a normalizable solution of (pp. 
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4.2.3 Reduction to large components. 

Finally, Breit used a procedure, which came to be known as "a reduction to 
large components" , to obtain from ( |16"D and (]XSf) the previously known forms 
of the operators which gave a 2 Ry level shifts. Recall first that in the stan- 
dard representation of the Dirac matrices it is convenient to organize the four 
components of the Dirac wave function of an electron into two Pauli 
spinors, ^u(l), the first and second (or upper) components of if), and ^(l), 
the third and fourth (or lower) components of For a hydrogenic bound 

state ipL is of order (v/c)ipu, so in this context the upper are the "large" 
and the lower are the "small" components of ^(l), respectively. One can 
similarly decompose a 16-component spinor if)(l, 2) into upper-upper compo- 
nents if)uu(l,2), lower-upper components i[)lu(1,2), etc. Breit wrote down 
the four coupled equations for the ipxY which follow from ([16]) and from three 
of these obtained approximations for ipm, ^Pul, and ipLL in terms of ipuu, 
the large-large components. When substituted in the remaining equation 
the result is an equation for ipuu which reduces to (|3]) in the extreme n.r. 
limit and contains further terms which can be identified as the previously 
known spin-self-orbit, kinetic energy, and spin- independent contact interac- 
tions. Further, on using the leading approximations to the small components 
(e.g. i/;lu(1,2)) ~ (ayp° p /2m)^ uu (l, 2)) in the expression (|18|), Breit found 
that AE reduced to the expectation value with the nonrelativistic wave func- 
tion of the orbit-orbit, spin-other-orbit, and spin-spin interactions mentioned 
above. 

4.3 Critique of the early work 

The most serious flaw of Breit 's work is the circumstance that the zero-order 
equation, the Dirac-Coulomb equation ((TH), has no bound states; that is, 
Hue has no normalizable eigenfunctions! This startling fact was first pointed 
out by Brown and Ravenhall, some 20 years after Breit's work.|| Its validity 
is most easily seen by first turning off Uc{r) in H^c- The resulting Hamilto- 
nian does have normalizable eigenfunctions, of the form ij)$ n = u m (l)u n (2), 
with u m (l) and u n (2) denoting normalizable eigenfunction of h e £ l {l) and 
/iff*(2), with enegies e m and e n , respectively. But such product states are 
degenerate in energy with similar product functions in which, say, "1" is in a 
continuum state with positive energy e\ > m and "2" in a continuum state 
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with negative energy e 2 < —m and e\ + e 2 = e m + e n . Since the sum e% + e 2 
can be kept constant with t\ increasing and €2 decreasing, the degeneracy 
is continuous. Thus, turning on any interaction which has non-zero matrix 
elements between positive- and negative-energy eigenfunctions of Hq Di such 
as e 2 /r 12 , will lead to a mixture of ipmn wftli a sea of non-normalizable states, 
i.e. to "continuum dissolution." |7]] 

The CD problem undermines both the derivation of the correction terms 
not associated with transverse photon exchange and the derivation of fll8|) , 
which purports to take these into account. However, suppose one overlooks 
the CD problem and i) considers the wave function ip(°\l, 2) with components 

^ 0) LU (l, 2) « (at ■ Pl /2m)4°i(l, 2), (20) 

etc., with ^^(1,2) identified with the n.r. SP wave function, as an approx- 
imation to an eigenfunction of an as yet unspecified relativistic zero-order 
Hamiltonian, (ii) accepts the formula (|T8| ) and evaluates it with the function 
ip(°\l,2), as was done by Breit. On carrying out some integration by parts 
and being careful to keep nonvanishing surface terms, one finds that, apart 
from the terms found by him, there is an additional term, a spin-spin contact 
interaction between electrons^, [J: 

= -(87r/3)(e 2 /47r) ( r 1 ■ a 2 5(r)m- 2 . (21) 

This term is simply the counterpart of the Fermi contact interaction between 
an electron and a proton, which had been known for many years previously in 
the context of the hydrogen hyperfine structure. It is somewhat scandalous 
that it took so many years to include this effect in fine-structure calculations 
in atomic physics. 

In summary, Breit 's theory of the helium atom was flawed both on con- 
ceptual and technical grounds. To a large extent the reason that this was 
not noticed for more than two decades is that he got the "expected answer" 
for terms previously known and that experiments were not precise enough to 
measure level shifts of order a 2 Ry with great accuracy. 

A sign that something was badly amiss with the general theory should 
have come from the following fact. Breit briefly considered the possibility 
of treating U B nonperturbatively by adding it to Uc and using as a starting 
point what I call the Dirac-Coulomb-Breit (DCB) Hamiltonian: 

H DCB (1, 2) = h%\l) + h e S\2) + U c {r u ) + U B (r 12 ) (22) 
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However, he soon noted that use of Ub in second-order perturbation theory 
led to strange results: The contribution from negative-energy states is readily 
seen be of order < U B /Am > and since U B ~ Uq this is of order a 2 Ry, as 
large as the first-order contribution of Ub*- So the Breit operator always came 
with the caveat: It was to be used only in first-order perturbation theory. 

Nevertheless, it should be emphasized that Breit's work was a pioneering 
example of the effort to derive effective potentials from quantum field theory, 
apparently accomplished when the first paper on the general quantum theory 
of fields by Heisenberg and Pauli was available only in preprint form. 

5 Quantum field theory of two-electron atoms: 
Foundations and level-shifts to order ce 3 Ry 



5.1 Foundations 

In the 50's, motivated by the increasingly accurate experimental results which 
were becoming available, e.g. for the splitting of the 2 3 P states and for the 
ionization energy, the foundation of the theory of the level-structure of He 



and ife-like ions was reexamined, independently by H. Araki|T0[ and me.|l 



The starting point for both these studies was a generalization of the two-body 



Bethe-Salpeter equation |T2| to include an external field, an integral equation 



for a 16-component spinor function ^(£1,22), of the form 

ip{x 1 ,x 2 ) = —i I K a (x 1 ,x 3 )K b (x 2 ,x 4 ;)G{x 3 ,x 4 ;x 5 ,x 6 )4>{x 5 ,x 6 )dT. (23) 



Here dr = d^x^d^x^x^d^x^ is the product of four-dimensional volume ele- 
ments, K a (xi,x 3 ) is the external-field Feynman propagator for electron "a", 
and G is an interaction kernel, constructed from two-particle-irreducible x- 
space Feynman graphs, involving the emission and absorption of one or more 
photons. In these graphs, the electron propagators are the K a and the pho- 
ton propagators depend on the choice of gauge. If one works in Coulomb 
gauge one can, by an extension of a method of Salpeter,|jl2[ derive an ex- 



act equation for the spatial part 0(xi,x 2 ) of the equal-times wave function 
ip(x.x,t; x 2 ,t) associated with a stationary state of energy E, 



^(x 1} t;x 2 ,t) =0(x 1 ,x 2 )e- iB '. (24) 
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This equation, which I shall not write down, is still linear but involves the 
eigenvalue E nonlinearly. On neglecting effects associated with electron- 
positron pairs and transverse photons in this equation, one arrives at what I 
call the "no-pair Coulomb- ladder (NPCL) equation" : 

H ++ (f)c(xi, x 2 ) = £0c(xi, x 2 ). (25) 

Here 

H ++ = ht{l) + fcS*(2) + L ++ U c (r 12 )L ++ (26) 
is the NPCL Hamiltonian and 



L++=L + (l)L + (2) (27) 



is the product of external- field positive-energy projection operators L + (i), 
defined by 

L+(i) = E | (j) n (i) >< <j) n {i) | . (28) 

Here the sum is over the positive-energy eigenstates <ft n (bound and contin- 
uum) of h% xt \i). The NPCL wave functions <fic of physical interest satisfy 
the constraint 

L + (i)<P = <P (z = l,2). (29) 
This is consistent with since the L + (i) commute with H ++ , 

[L + (i),H ++ ]=0. (30) 

The nomenclature "NPCL" arises from the fact that (^) effectively sums 
all time-ordered diagrams which involve repeated Coulomb scattering of the 
electrons in the presence of the external field, with no electron-positron pairs 
present in intermediate states. 



The NPCL equation fl25|) has the following features: (a) Unlike the DC 
equation (^) it does not suffer from continuum dissolution, because when 
Uci^u) is turned on the projection operator L ++ prevents the mixing referred 
to above, (b) It is exact to all orders in the external field strength, (c) 
Because of the condition (p9|) it can be reduced without approximation to 
SP form, with the eigenvalue E still appearing linearly, (d) This relativistic 
SP equation reduces to (0) in the extreme NR limit and contains all the 
a 2 Ry corrections not associated with transverse photons. 
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The NPCL wavefunction can therefore be used as the basis of a pertur- 
bative treatment of the four- dimensional equation (|23|) , as developed in Ref. 
11, to compute effects associated with transverse-photon exchange and vir- 
tual pairs. In particular, the leading effects of transverse-photon exchange, 
of order a 2 Ry, can be incorporated by simply adding Ub to Uc in ([261) : 

H' ++ = h%\l) + fcS*(2) + L ++ [U c (r 12 ) + U B (r 12 )]L ++ ). (31) 

The corresponding equation 

H' ++( j>' c = E'(j)' c (32) 

has, like (|25|), a clear-cut origin in field theory: It is obtained from ( p3|) by 
keeping the one-photon exchange graph in the interaction kernel G, but ne- 
glecting the time-component ko in the factor (k 2 + ie)^ 1 appearing in the 
transverse-photon propagator. It can be regarded as a QFT-based replace- 
ment for the patchwork approach to fine structure. The eigenvalues E' in- 
corporate all fine-structure effects of order a 2 Ry and of course some of the 
effects of order a 3 Ry and higher. 

The basic difference between the modern approach and the earlier one is 
not the four-dimensional formalism, which is mainly a convenient tool and 
can to a large extent be dispensed with, but that QED incorporates the 
ideas of Dirac's hole theory rather than the Dirac's one-electron theory. The 
appearance of the projection operators is the mathematical expression, in 
the present context, of the idea that the negative-energy states are filled and 
transitions to them are forbidden by the Pauli exclusion principle. 

5.2 a 2 Ry and a 3 Ry level-shifts 

Analysis of (|26|) shows that the low-lying eigenvalues have the form 

E = 2m + W + 5W C (33) 
where W is the n.r. eigenvalue of (0) and 

SW C = 5W { c ] + 5W [ c ] + ... (34) 
The first term in ([34]), of order a 2 Ry, is defined by 
5W [ c ] =< (j) | -"pi/8m 3 " + (Tta/2m 2 )[Z5(r 1 ) - 5(r 12 )} + 1 <-> 2 | > (35) 
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+ < I (a/4m 2 )ai • (Zrx x Pl /r 3 - r 12 x pi/r? 2 ) + 1 «-> 2 | > (36) 

and is the sum, for each electron, of a kinetic energy correction, a spin- 
independent contact interaction, and a spin-self-orbit interaction. The latter 
terms are natural extensions of those found from the Dirac equation in the 
one-electron case, if one recognizes that electron "2" serves as the source of 
an electric field just as well as the nucleus. The quotes on the p 4 term serve 
SIS db reminder that this term must be evaluated as 

- < P?0 I P?0 > /8m 3 , (37) 
to avoid double counting. [|] The second term in (|34"D, defined by 

5W§ ] = -(a 2 /m 2 )(ir/2 + 5/3) < | <5(r 12 ) | >, (38) 
is an s.i. contact term of order a 3 Ry, arising directly from 



The remaining corrections of order a 2 Ry and another correction of order 
a 3 Ry may be obtained by a similar analysis starting from (|3l"D or, equiva- 
lently, by evaluating 

5E B =< <p c \U B \<j)c> (39) 
with a sufficiently accurate approximation to 0c- The result is 

5E B = 5W B 2) + 5W { b ] + ... (40) 

The first term in (|40|) which comes from reduction of the Breit operator, 
using approximations such as (|2"0"D, is given by 

5W B 2) =< | C/ _ + U m + U s „ s | > (41) 

with 

Uo- = (a/2m 2 )rr 2 1 [ P rp 2 + r-(f ■ Pl )-p 2 ] (42) 
a form of the orbit-orbit interaction, 

U s - - = -(a/2m 2 )r u 3 [a r r 12 xp 2 +l«2] (43) 

the spin-other-orbit interaction, and U s - S the full spin-spin interaction, in- 
cluding the contact term (21): 

E/,_ a = (a/Am 2 )[r^(a v a 2 - 3a r fa 2 -f) - (tor/3)<T V <T 2 6(r)]. (44) 
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The second term in ( fiOD is a spin-dependent partner of (|38| ) and contributes 
to fine-structure splitting in order a 3 Ry: 

Swf = -(a 2 /^ 2 )(4/3)(vr/2 + l)<4>\\a 1 - a 2 6(r 12 ) \<f>>. (45) 



The corrections (|38|) and (f45|) are both associated with the fact that the usual 
approximations to the LU, UL and LL components of the 16-component wave 
function are not sufficiently accurate at short distances, for the purpose at 
hand. 

The computation of the remaining terms of order a 3 Ry requires a return 
to (|23|) and i) keeping k in the transverse photon propagator (this corre- 
sponds to taking into account the Coulomb interaction during the exchange 
of one transverse photon), ii) keeping the irreducible two-photon exchange 
graph in the kernel G, and iii) keeping graphs associated with radiative cor- 
rections, the latter of course include self-energy effects, associated in hydro- 
gen with the Lamb-shift. In He, such effects were first estimated by Kabir 



and Salpeter.[13|] So much for the 50's. Later work in this area will be men- 
tioned in Sec. 8. 



6 Radiative decays of helium: Enter the Sun 
once more 

Because spin-dependent forces are weak if Z is not too large, the states of 
He or low-Z ife-like ions may be classified into singlet and triplet states, 
corresponding to total spin and 1, respectively. The lowest lying triplet 
S-state, described as 2 3 S' 1 , lies above the l 1 ^ ground state and is expected 
to be very long-lived, because of the form of the (nonrelativistic) radiation 
operator for one-photon decay, viz. 

Rrad = (e/v / 2cJ)S(p ii /m — isi x k/m) • e*exp(— ik ■ r-j). (46) 

The orbital part of R ra d vanishes between L = states and in the dipole 
approximation the spin part of R ra d becomes proportional to the total spin S. 
So the matrix element for the would-be magnetic dipole transition (A J = ±1, 
0; no parity change): 

2 3 5 1 ^l 1 S + 7 (47) 
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is doubly vanishing in this approximation: The matrix element of S is zero 
between a triplet and a singlet spin wave function and the spatial parts 
of the initial and final states, belonging to different eigenvalues of (3), are 
orthogonal. 

What then is the rate for this decay? This problem was apparently first 
studied in 1940 by Breit and Teller (BT) who, motivated by questions in 
astrophysics, wrote a pioneering paper on the I7 and 27 decays of H and 
He.^Q They estimated that R(2 3 S 1 S + r) « 10~ 24 i? where R denotes 
the rate expected for an El decay of unit oscillator strength (« 10 9 s~ 1 ), 
corresponding to a lifetime for the I7 decay of 10 8 years! This was the last 
word on the subject for almost three decades, when the Sun enters our story 
once more. In 1969 it was noted by A. H. Gabriel and C. Jordan|I3J that 
certain lines seen in the soft x-ray spectrum of the solar corona correspond 
in wavelength to the transition ( pT7| ) in certain He-like ions, ranging from C 
V to Si XIII, and suggested that these transitions were in fact responsible for 
the observed lines. However, they noted that the intensities for these lines 
would be extremely small, if the result of BT were simply scaled up, even 
with a high power of Z. 

Soon thereafter H. Griem pointed out that BT vastly underestimated 
the rate for (|47|) in He and made an estimate, patterned on the (largely) 
correct calculation of BT for the hydrogenic decay 2 1//2 S , i/ 2 — > 2 1 / 2 Si/ 2 + I7, 
which indicated that the Gabriel- Jordan interpretation was plausible. [|16j A 
similar estimate was made by G. Feinberg. |L7] At about the same time the 
whole subject was brought down to Earth, so to speak, with the development 
of beam foil spectroscopy and the terrestrial measurement of such rates, 
beginning with that of (fi? ) in Ar 16+ by Marrus and Schmieder.||18|| This 
raised interest in the whole subject and led to the development of a general 



QED based theory of such decays. |L9L 20 



A correct starting point for the calculation of one-photon decays in if -like 
ions which includes relativistic effects is the equation 

Mi =< -0/ I a ■ e*e" ikr | fa > . (48) 

where the ip 's are normalized solutions of (p~0|) . This is also what one would 
write down in Dirac's one-electron theory and it yields the correct leading 
terms for all such decays. Thus the counterpart of ( pE8| ) for a two-electron 
system is 

M 2 = £ < V/(l, 2) I a r e*e-^ ri | fa(l, 2) >, (49) 
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where now the ^'s are solutions of (p5|), which is according to the above 
discussion a QED-based analogue of ([48|) . However, it turns out that (|4^) 
does not give the full leading terms for the Ml decays in question: One 
must go back to QED and include time-ordered graphs in which, e.g., the 
Coulomb field of electron "1" or the nucleus creates a virtual e~ — e + pair 
and the positron annihilates with electron "2" to create the emitted photon. 
So field theory rears its head again. The final result is 

M(2 3 S 1 -> l'So + 7) ~ 2zS /4 T /4 (50) 

where 

T fi =< <f> f I p 2 /3m 2 + k 2 r\/l2 - r x U\rx)/&m \ fa > (51) 

and 

Z fi =<Xf\o-i-kxe*/m\xi> (52) 

with Xi and Xf the triplet and singlet spin wave functions. The corresponding 
decay rate is then found to be 

Rue-uke = (2/3)(8ak 3 /m 2 ) | T fi | 2 (53) 

The evaluation of Tfi for large Z was first carried out by Drake using many- 
parameter variational wave functions and an expansion in powers of 1/Z. 
Because of early discrepancies in the cases of Ar l6+ and C7 15+ , further work to 
estimate and include radiative, recoil, retardation and relativistic corrections 



were carried out by various workers. [|2T| It turns out that for Z not too small, 



larger than say 15, these results can be put together to write a cooperative 
explicit formula: plj 



Rue-iike = (2/3)«(aZ) 10 (m/972)[l - A.10/Z + 6.7 /Z 2 + l.07(aZ) 2 ] (54) 

which could, conservatively, be expected to be accurate to better than 2% 
in the range 15 < Z < 40. Comparison with data ranging from sulfur (Z 
= 16) to krypton (Z = 36), for which the lifetime ranges from 700 ns to 
.2 ns, a factor of more than 3000, indeed showed agreement of theory and 
experiment to within a few percent in all cases. When one includes He, for 



which the lifetime of the 2 Si state is about 10 s, 22 albeit with a rather 



large (30%) error, and evaluates Tfi with a Hylleraas type of wave function, 
the agreement extends over a range covering 14 orders of magnitude. This 
was surely a triumph of QED in a rather unusual domain. Again, recent 
work in this area will be noted in Sec. 8. 
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7 Retarded van der Waals forces and Ryd- 
berg states of helium. 

Van der Waals forces are, in the first instance, the (spin-independent) elec- 
tromagnetic forces which act between atoms at separations R large compared 
to atomic size a. Two years after the development of quantum mechanics, 
S.C. Wang and later F. London showed that the electrostatic interaction 
between the constituents of two neutral atoms gives rise to an effective inter- 
atomic potential of the form U = —C / R 6 for R ^> a (often called the London 
potential).^] But it took two decades after the development of QED, for 
H.B.G. Casimir and D. Polder to point out that the inclusion of the effects 
of transverse photon exchange between the electrons leads to a potential 
which falls of as R~ 7 at very large separations, of the order of the maximum 



wavelength for electric dipole excitation. |24j Inclusion of transverse photon 
exchange corresponds, in the framework of QFT, to taking into account the 
retarded character of electromagnetic interactions in CED. Thus one refers 
to the potential at very large R either as the Casimir-Polder potential or as 
the retarded van der Waals potential. Unfortunately, because the potential 
is very small in this region, direct experimental detection of the retarded 
character of atom-atom van der Waals forces has proved to be elusive. 

On general grounds, one expects that related retardation effects manifest 
themselves between a polarizable neutral system A and a charged system 
B, or between two charged systems A and B. However in these cases, such 
effects will only be corrections to the dominant forces, corresponding to the 
electric dipole potential Udi P = — a^e^/2i? 4 and the Coulomb potential Uc, 
respectively. It was pointed out by L. Spruch that precise measurement of the 
fine structure of Rydberg levels in helium might be able to detect such cor- 
rections and S. Lundgren and co-workers undertook a series of experiments 



to detect such effects, specifically on the n — 10 levels. |25[ This stimulated a 
great deal of further work on the QED aspects by theorists including Spruch 
and J.F. Babb, and C.K. Au, G. Feinberg and myself. [p5|] There is only one 
point I will elaborate on in this regard, by way of making a connection with 
the beginning of the modern theory of helium. G.W.F. Drake, by using a 
doubled set of Hylleraas-type wave functions and modern computer power, 
has extended the variational methods pioneered by Hylleraas so long ago to 
achieve unprecedented accuracy in computing the eigenvalues and eigenfunc- 
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tions of the nonrelativistic Hamiltonian (3).|[25|| R. Drachman has achieved 
very high accuracy for the large n, large / states, by using a systematic per- 
turbation expansion in the residual interaction e 2 (l/r 12 — 1/Vi) of the outer 
electron with the core.|25[ Thus helium now approaches hydrogen in hav- 
ing an exactly-solvable zero-order problem. As a result, the calculation of 
relatively high-order relativistic and QED type corrections in helium and 
helium-like ions has become much more meaningful than in earlier decades, 
as emphasized below. 



8 Some recent developments; Concluding re- 
marks 

As described in Sec. 5, the order a 3 Ry corrections to the low-lying a 2 Ry fine 
structure of helium were calculated long ago. Increased experimental preci- 
sion and the possibility of using He fine-structure for obtaining values of the 
fine-structure constant competitive with other methods, has led to improve- 
ments on the earlier calculations and to the computation of still higher-order 
effects. Helium has also been used as a testing ground for techniques vital for 
studying atoms with more than two electrons, such as relativistic many-body 
perturbation theory (MBPT), configuration interaction (CI) methods, and 
multi-configuration Dirac-Fock (MCDF) methods. 

A major step in the direction of going to higher orders was taken by M.H. 
Douglas and N. Kroll, who, using the methods of Ref. 11, obtained formulas 
for fine-structure splitting of order a 4 i??/.|26| Recently this work has been 



improved and extended by several theoretical groups. T. Zhang and G.W.F. 
Drake pTfl and Eides et al. have verified the correctness of the somewhat 



phenomenological treatment of radiative corrections by Douglas and Kroll. 
Z.-C. Yan and Drake [^] have computed the fine-structure of the 2 3 P state 
in He and low-Z He-like ions to a few parts in 10 9 , obtaining good agreement 
with recent high precision measurements PUJ for Z = 3,5, and 9 but not for 



Z = 2 and 4. Complementary calculations for higher Z have been carried out 
by M.H. Chen et aZ.,|3l] using relativistic CI methods, and by D.R. Plante et 



a/.|32[, using relativistic MBPT methods. In a heroic computational effort, 
T. Zhang has obtained corrections of order aHnaRy to energy levels and 



corrections of order a InaRy and a Ry to fine- structure splitting. [ 33| These 
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will be important in using helium fine structure to obtain an independent 
determination of a; a recent measurement by Myers et al. [34] of the fine- 
structure splitting of the 2 3 P splitting for Z = 7 is of sufficient accuracy to 
be sensitive to such corrections. 

Much progress has been made in recent years in developing new tech- 
niques for the accurate calculation of self-energy corrections in He and He— like 
ions, i.e. corrections which are the counterpart of what is referred to as the 
Lamb shift in if and if -like ions. A recent example, where further references 
may be found, is the work of Persson et al. [p5|| , in which good agreement is 
found with experimental results of Marrs et a£.[56 on several iie-like ions, 
ranging from Z = 32 to Z = 83. 

With regard to decay rates, W.R. Johnson et aZ.|3"7|l have carried out 
extensive relativistic calculations of radiative decay rates in the helium iso- 
electronic sequence, using the no-pair Hamiltonian (26) as a starting point. 
These calculations are so accurate that disagreement with experimental val- 
ues could be interpreted as associated with radiative corrections to decay 
rates, a relatively unexplored subject. Very recently, P. Indelicato|!8| has 
extended the MCDF method to the calculation of the forbidden magnetic- 
dipole decays, discussed in Sec. 6. 

As a final example, in a quite different direction, of the way in which 
helium engages the attention of physicists, I refer you to a paper by Deilamian 



et al. ||39|| , who have made a search for He states which are forbidden by 
the Pauli antisymmetry principle, by way of putting an upper bound on its 
hypothetical violation. 

The story goes that when a mathematics professor was asked by a student 
"What is the use of number theory?", he replied: "Number theory is useful 
because you can get a Ph.D. with it." As one may infer from the above saga, 
helium has certainly been useful for getting a doctoral degree, and its uses 
in the laboratory and in providing impressive views of football stadiums are 
well known. In a more serious vein, helium, together with its ionic partners, 
is likely to continue for a long time to come its unique role as a laboratory 
for the application of quantum field theory and the techniques of mathe- 
matical physics to a system of fundamental physical interest. It provides a 
unique setting for the study of the electron-electron interaction, the electron 
being one of the few observed particles which are still usefully regarded as 
elementary. While helium may be inert chemically, it remains highly active 
intellectually. 
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